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Abstract Combining two-color infared pump-probe spectroscopy and anharmonic 
force field calculations we characterize the anharmonic coupling patterns between 
fingerprint modes and the hydrogen-bonded symmetric Vnh2 stretching vibration 
in adenine-thymine dA2o-dT2o DNA oligomers. Specifically, it is shown that the 
anharmonic coupling between the 5nh2 bending and the Vc4=04 stretching vibration, 
both absorbing around 1665 cm~^, can be used to assign the Vnht fundamental 
transition at 3215 cm^' despite the broad background absorption of water. 



1 Introduction 

Vibrational energy redistribution and relaxation in complex systems depends on the 
network of anharmonically coupled vibrational states subject to fluctuations due to 
the interaction with some environment [ 1 1. Focussing on hydrogen-bonded systems 
there is considerable evidence that the time scales for relaxation can be in the subpi- 
cosecond range pointing to a rather strong interaction, e.g., of the excited stretching 
vibration with other hydrogen bond (HB) related modes such as the bending and the 
low-frequency HB distance vibration as well as with the solvent |l2][3]|4]|5]|6l- 
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One of the most prominent hydrogen-bonded systems is DNA. Despite numerous 
experimental and theoretical investigations on vibrational spectra of nucleic acid 
bases ||7]|8]|9][T0][TT1[12][T3], information on inter- and intramolecular interactions 
in base pairs and DNA oligomers is still limited ll4l[T5l[T6l[T7l[T8l[T9ll20ll2n 
I22I I23I |24] |25l. A recent example is the work on single adenine-uracil (AU) base 
pairs in the Watson-Crick geometry in solution, which showed an enhancement of 
vibrational energy relaxation of the NH stretching vibration by a factor of three as 
compared to the isolated uracil base 1211 . 

DNA oligomers adopt different types of conformations, both in gas and con- 
densed phases, such as the A, B, B', C, D, and Z form, depending on water and salt 
concentration, type of cations, pH, and base sequences ll8l [T2ll25ll26ll27ll28ll29ll . In 
the condensed phase the conformations of DNA oligomers are stabilized by wa- 
ter molecules that form water networks, predominantly in the major and minor 
grooves, and near the phosphate groups of the backbone IIT2I . Among the differ- 
ent types of base sequences, adenine-thymine (AT) oligomers are special because 
they do not undergo transitions from the B to the A form upon reducing the water 
content. Instead, AT oligomers adopt the B' form at low water concentrations, with 
4 to 6 water molecules per base pair that may be hydrogen-bonded to the oligomer 
inilSOlllllllll. In the B' form of the AT DNA oligomer two HBs are formed in the 
Watson-Crick configuration, i.e., between oxygen (04) of the thymine and the NH2 
group of the adenine (N6), and between the NH group of the thymine (N3T) and the 
nitrogen atom of the adenine (Nl), see Fig.[T] 

Vibrational modes expected to be strongly influenced by the hydrogen-bonding 
in the DNA helix are the carbonyl stretches Vc2=02at 1716 cm^' and Vc4=o4at 
1665 cm-i and the amine bending 5nh2 at 1665 cm-^ E] [H H] [HlMl |35l |36l 
[37I l38l . Note, that in contrast to H2O, in D2O the 5nd9 vibration of adenine and 
the carbonyl vibrations of thymine are decoupled, due to the frequency shift from 
5nh2 to 5nDt II22I |23l . The 5htO vibration of water molecules in DNA samples 
typically absorbs in the same spectral region, i.e., around 1650 cm~' ||7l |35] [39l . 
A direct experimental assignment of Vnht and Vnh in AT DNA oligomers in the 
condensed phase is very difficult. Typically, symmetric and antisymmetric Vnh2 
stretching vibrations absorb around 3300 cm^' 121. However, the spectral range 
from 3050 to 3600 cm^ ' is dominated by the strong absorption of the water OH 
stretching vibration. Reducing the water content of the DNA oligomers does not 
solve this problem, because at extremely low water contents the DNA oligomers do 
not adopt a well defined structure. 

Ultrafast time-resolved infrared (IR) spectroscopy is ideally suited to address this 
issue as has been shown in studies of inter- and intramolecular couplings and energy 
relaxation dynamics in various hydrogen-bonded systems 121 [3l |4l |40l l4ll . In this 
contribution we focus on shifts in oligomer vibrational modes induced by excitation 
of the Vc2=02 or the Vc4=o4 / ^nht oligomer fingerprint vibration. These shifts orig- 
inate from inter- and intramolecular couplings among different vibrational modes of 
the DNA oligomer and depend on the strength of the couplings as well as the energy 
mismatch between different transitions. Related effects are particularly pronounced 
if overtones or combination modes match a fundamental vibrational transition (res- 




Fig. 1 (a) Scheme of the AT DNA Watson-Crick configuration, b) Structure of a single AT DNA 
base pair in the Watson-Crick configuration of a dodecamer (taken from 428d.pdb 1331 ). The oxy- 
gen atoms of water molecules forming HBs to the nucleic acids are presented as spheres. Distances 
of the oxygen atoms of water molecules in the major groove a, b, and c to the 04, N6, and N7 atoms 
are 2.93 A, 2.93 A, and 2.80 A, respectively. Water molecules d, e, and f of the minor groove have 
distances of 2.88 A, 4. 16 A, and 2.61 A to the 02 and N3 atoms, respectively. Intrastrand distances 
of the 02, 04, and N6 atoms to neighbooring thymine and adenine bases are 4.12 A, 3.53 A, and 
3.34 A, respectively. 



onance enhancement). This already affects the linear absorption band shape, but 
also the vibrational relaxation dynamics [jSj. A particular strength of the ultrafast IR 
pump-probe spectroscopy is the capabiUty of uncovering vibrational spectral fea- 
tures not visible in linear spectroscopy due to excessive solvent absorption. This is 
demonstrated in the experiments presented here, where we excite oligomer vibra- 
tions between 1600 and 1760 cm^' and probe for the oligomer Vnh2 vibration in 
the region of 3050 - 3250 cm^ ' , which is dominated by water absorption. 

The experimental assignment of the adenine Vnh2 vibration and the coupling 
pattern across the HBs is supported by quantum chemical calculations of anhar- 
monic couplings which are used for obtaining fundamental transitions frequencies 
for a set of relevant modes of a microsolvated gas phase AT model. In principle 
accurate theoretical modelUng of the vibrational dynamics of DNA AT base pairs 
requires taking into account several effects: (i) The intermolecular double HB be- 
tween adenine and thymine, (ii) The interaction between different base pairs along 
the DNA strand, (iii) The charges as well as the dynamics of the backbone, and (iv) 
The influence of water molecules which may, for instance, make a HB to the base 
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pair. Here, we are aiming to obtain a semiquantitative understanding of the transient 
band shifts, whereby it is assumed that they are dominated by effect (i), that is, the 
anharmonic coupling pattern due to the intermolecular HB. The effect of (ii)-(iv) 
can be of static nature, e.g., changes in the anharmonic frequencies and coupling 
constants, and also of dynamic nature, e.g., fluctuation of the energy levels. How- 
ever, here we will only focus on the static influence of a well-defined environment 
determined by microsolvation of the AT base pair by several water molecules. 

Isolated and microsolvated base pairs have been extensively studied theoretically, 
focussing in particular on the stability of different isomers, see, e.g. the work by 
Hobza and coworkers 1421 l43l l44l as well as by Fonseca-Guerra et al. B31 . Al- 
though there is a number of reports on potential energy surfaces of base pairs in 
harmonic approximation, there appear to be only a few calculations addressing an- 
harmonicity in the context of, e.g., proton transfer ifTol l46l l47l . the coupling to 
the intermolecular HB vibration |15| or the assignment of different gas phase iso- 
mers 1481 . Most notable in this respect is the recent study of the anharmonic spec- 
trum of a guanine-cytosine pair P9|| as well as the development of a vibrational 
exciton model to describe nonlinear IR spectra involving DNA fingerprint modes 

This Chapter is organized as follows. In the next Section we will first discuss 
the effect of solvating water molecules on the anharmonic IR spectra of an isolated 
AT pair. In this context we will scrutinize the applicability of a dual level approach 
which combines different quantum chemistry methods within a correlation expan- 
sion of the potential energy surface (PES). For the case of two water molecules we 
will present an analysis of the anharmonic coupling patterns between the Vc2=02 , 
Vc4=04and 5nh9 vibrations and the symmetric Vnht mode in Section l272l Section 3 
gives details on the experimental setup and presents results of two-color pump-probe 
spectra. Finally, we give a comparison between theory and experiment in Section|4] 
which leads us to the assignment of the Vnh2 fundamental transition. 



2 Microsolvated AT Base Pairs 

2.1 Fundamental Transitions Using a Dual Level Approach 

In the following we present results on fundamental vibrational transitions of isolated 
AT base pairs microsolvated with 1-4 water molecules. The aim of this study is 
twofold: First to find out about overall changes of IR transitions of base pair modes 
due to the interaction with water molecules. And, second, to test the performance 
of a dual level approach combining density functional (DFT) and semiempirical 
PM3 data to expand the PES. Throughout we will assume that the deviations from 
equilibrium structures are small enough such as to allow the use of normal mode 
coordinates Q for spanning the PES, i.e. V = V{Q). 
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Table 1 Harmonic frequencies (in cm^' ) for the target modes of the different model structures 
as obtained using DFT/B3LYP with a 6-31++G(d,p) basis set. The displacement vectors for the 
solvated structures are shown in Figs.|2]and|3] 



modes 


AT 


AT-H2O AT-(H20)2 


AT-(H20)3 


AT-(H20)4 


Vnh, 


3410 


3393 


3401 


3415 


3410 


VC2=02 


1797 


1794 


1799 


1799 


1776 


5nH2 


1689 


1718 


1720 


1727 


1727 


VC4=04 


1728 


1731 


1714 


1720 


1720 



Under the conditions of the experiment there are 4 to 6 water molecules per AT 
pair which can form different HBs to the base pair (see also Fig. [U. Our interest 
will be in the IR transitions of the NH2 and C4=04 groups such that water situated 
in the major groove shall be of importance. However, for comparison we also con- 
sider a structure where a water molecule is on the C2=02 side. There are several 
microsolvation studies which focussed on the effect of water on base pair properties 
such as interaction energies or HB lengths 1431 l44l |45l |5T| . To our knowledge there 
is, however, no theoretical account on anharmonic IR spectra of HB related modes. 
The four structures which will be discussed in the following are shown in Figs. |2] 
and[3] They have been obtained by geometry optimization at the DFT/B3LYP with 
a 6-31++G(d,p) level of theory using Gaussian 03 1521 . Notice that these are not 
necessarily the most stable structures at this level of theory (see also discussion in 
Ref. El). Our choice has been biased by the requirement that the water molecules 
should be close or even hydrogen-bonded to the considered target modes. The latter 
are shown in terms of their normal mode displacement vectors in Figs. |2] and [3] as 
well. The respective harmonic frequencies are compiled in Table [T] 

In AT-H2O, Fig.|2](left column), the water molecule is hydrogen-bonded between 
the adenine N6-H and the N7 sites. This causes the Vnh2 vibration to acquire some 
water stretching character lowering its harmonic frequency. The 5nh, vibration is 
only slightly mixed with some water motion and essentially constrained so that its 
frequency is blue-shifted. The next water molecule in AT-(H20)2, Fig.|2](right col- 
umn), makes a HB to the oxygen of C4=04 lowering the Vc4=o4 frequency slightly, 
but at the same time mixing this vibration with 5nh9 type motions. For the case of 
three water molecules, Fig.|3](left column), there is the possibility to form a hydro- 
gen bonded water chain connecting the 04, N6-H, and N7 sites. This reduces the 
mixing of the Vnh2 ™d water motions, but the 5nh, vibration contains a water bend- 
ing now as does the Vc4=o4 mode. Adding another water at the C2=02 site leads as 
expected to a shift of the Vc2=02 transition only. Fig. [3] (right column). Overall we 
notice that the presence of solvating water molecules has the strongest impact on the 
5nh, and Vc4=04 vibrations, with the latter acquiring substantial 5nh2 character 

So far we have only discussed harmonic frequencies. The effect of anharmonicity 
can be treated using either a Taylor expansion of the PES in terms of normal mode 
coordinates or by explicitly spanning the PES on a numerical grid. The discussion 
of anharmonic force constants is postponed to the following section. Here, we will 
focus on an explicit PES generated by means of the following correlation expansion, 
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v{C4=04) 




Fig. 2 Normal mode displacement vectors for the target modes as calculated in harmonic approx- 
imation of the AT-(H20)„=i,2 PES using DFT/B3LYP with a 6-31++G(d,p) basis set. 



here written up to three-mode correlations, 1531 

I i<j i< j<k 

Neglecting rotation ll54l . i.e. assuming that the kinetic energy operator is diagonal, 
the eigenstates of the respective Hamiltonian can be obtained by straightforward di- 
agonalization using, e.g., the Lanczos method. For systems of the size of solvated 
base pairs the calculation of higher order correlation terms becomes rather expen- 
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Fig. 3 Normal mode displacement vectors for the target modes as calculated in harmonic approx- 
imation of the AT-(H20)„=3,4 PES using DFT/BSLYP with a 6-31++G(d,p) basis set. 



sive. Here, an interesting alternative are so-called dual level schemes where low- 
order correlation PES are calculated at a higher level of quantum chemistry than 
multi-mode correlation PES. For instance, Scheurer and coworker have combined 
MP2 and PM3 calculations to find a rather good description of IR spectra of a model 
peptide ll55l . 

In Table 12] we present results of dual level calculations on 4D models including 
the target modes of Figs. |2]and|3] Here the one-mode potentials, V''^(2;), have 
been calculated using the DFT method, while two- and three mode PES were gener- 
ated using the semiempirical PM3 approach. Comparing these anharmonic results 
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Table 2 Anharmonic frequencies (in cm^ ' ) for the target modes of the different model structures 
as obtained using Eq. (T} with the one mode potential generated by the DFT/B3LYP method with 
a 6-31++G(d,p) basis set and the two- and three-mode PES obtained by the PM3 approach. 



mode AT-H2O AT-(H20)2 AT-(H20)3 AT-(H20)4 



VnHt 


3326 


3332 


3307 


3310 


VC2=02 


1803 


1813 


1797 


1785 


5nH2 


1732 


1760 


1777 


1753 


VC4=04 


1701 


1664 


1635 


1643 



Table 3 Comparison of anharmonic frequencies (in cm ' ) for the target modes of the AT-(H20)2 
model structures of different two-dimensional calculations and different levels of quantum chem- 
istry (DFT: full DFT/B3LYP, 6-31-H-G(d,p); PM3: full PM3; DUAL: DFT one-mode and PM3 
two-mode PES). 



2D model 


mode 


DFT PM3 DUAL 


Vnh2 Ah, 


VnH2 

Snh-, 


3311 3154 
1727 1950 


3321 
1730 


VC2=02 , 5nH2 


VC2=02 


1796 1901 
1742 1703 


1796 
1742 


VC4=04 >5nH2 


VC4=04 
5nHt 


1715 1807 
1744 1671 


1680 
1777 



with the harmonic values in Table[T]we notice the following: The Vnh, vibration is 
strongly affected and red-shifts by about 70- 1 00 cm^ ' depending on the cluster size. 
The Vc2=02 vibration is only slightly affected. Essentially, these two modes behave 
as expected. Except for the AT-H2O case the 5nh2 vibration is blue-shifted by about 
40 cm^ ' . The fundamental transition of the Vc4=04 mode, on the other hand, is 
considerably red-shifted by about 30-80 cm^^ . In terms of the experimental assign- 
ment given in Table 6 below the Vc4=o4 agrees fairly well with experiment which 
puts this transition at 1665 cm^' Il7][8l[l4l. However, the 5nht vibration is believed 
to absorb around 1665 cm^' as well, which is at variance with the prediction of the 
dual level scheme. Moreover, in the harmonic case but also in the fourth order anhar- 
monic force field calculations reported below the frequencies of Vc4=04 and 5nh2 
are almost identical. It would be surprising if the higher order anharmonic terms 
implicitly included in the PES expansion changes this situation to such an extent. 

Since multimode calculations are rather expensive we have chosen to scrutinize 
the effect of the PM3 approximation by comparing different 2D models of AT- 
(H20)2 at the full DFT, PM3, and dual level. The results are compiled in Table 
[3] The correlation between the Vnht and 5nh2 modes is rather well-desribed by 
the dual level scheme, the error being just a few cm^ ' . The same holds true for 
the correlation between the Vc2=02 and 5nht modes. However, the hybrid scheme 
is performing poorly for the correlation between the Vc4=o4 and 5nh, modes and 
gives the frequency shifts of opposite sign also observed for the full 4D calcula- 
tion in Table |2] The failure of the dual level scheme to describe the coupling of 
the Vc4=04 and 5nh2 bending motion can be understood in terms of their consider- 
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able mixing as quantified by a normal mode internal coordinate decomposition. In 
general, an internal coordinate may contribute to the displacement along several nor- 
mal modes. Analyzing the present situation using the scheme of Boatz and Gordon 
ll56ll the internal coordinate describing the NH2 angle contributes to the decompo- 
sition of the 5nh2 bending normal mode by only 29 %, while its contribution to the 
Vc4=04 stretching normal mode (see also Fig.|2]) and several purine ring deforma- 
tions normal modes ranges between 15 and 17% . On the other hand, the internal 
NH2 angle does not contribute to the decomposition of the Vc2=02 and Vnht normal 
modes which are 92% stretching of the C2=02 bond and 82% and 17% stretchings 
of the N-H bonds, respectively. 



2.2 Anharmonic Coupling Patterns 

In this section we explore the second possibility to generate multidimensional PES, 
i.e. a Taylor expansion in terms of normal mode coordinate with respect to the ge- 
ometry of the stable structure. Including terms of to fourth order we have (using 
dimensionless coordinates) 

^- i ijk ^- ijkl 

Third and fourth order anharmonic coupUng constants are then calculated using a 
combination of analytical second derivatives and finite differences ll57l . Specifically, 
we have used the symmetric expressions 



^(3) _ +)iKij -iiK-j +Kij 

K^Z = — ■ (4) 

Here, Kf- is the Hessian calculated at displaced geometries where for the displace- 
ment we used AQ ^ 0.03 for the cubic and AQ~ 0.04 for the quartic force field. 
Note that for the construction of the Hamiltonian we have neglected contributions 
which are off-resonant by more than ~1000 cm^' . 

In the following we will focus on the case of two water molecules only, i.e. 
AT-(H20)2, since this already contains the essential effect of hydrogen-bonding wa- 
ters as discussed in the previous section. The diagonal force constants for the four 
target modes of Fig. |2] (right panel) are given in Table 4. Important third order an- 
harmonic coupling constants involving the Vnh9 mode are compiled in Table |5] As 
expected the diagonal anharmonic force constants are largest for the Vnht mode. 
More interesting, however, is the coupling pattern between this mode and the fin- 
gerprint modes. Here, we observe two dominating Fermi-type resonance couplings; 
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Table 4 Harmonic and diagonal anharmonic force constants (in cm ' ) of the relevant system 
modes. 





VnHt 


VC2=02 




Vc4=04 


CO 


3401 


1799 


1713 


1720 




-261 


61 


-19 


-31 




24 


5 


4 


2 



Table 5 Third order coupling constants between the Vnh2 mode and relevant fingerprint modes 
(in cm^' ) in AT-(H20)2. The numbers in parentheses refer to the isolated AT case. 





VC2=02 


5nHt 


VC4=04 


VC2=02 


-1 (-2) 


5(8) 


0(-2) 






96 (87) 


-66 (-17) 


VC4=04 






11 (2) 



Table 6 Anharmonic frequencies (in cm ^ ) for the target modes of the AT-(H20)2 model calcu- 
lated usino the anharmonic expansion, Eq. (|2j and the DFT/B3LYP method with a 6-31-l~l-G(d,p) 
basis set. The 4D model is compared with a 6D model which includes in addition the most strongly 
coupled water stretching and bending modes at the N6-H site whose anharmonic frequencies are 
3752 cm^' and 1588 cm^' , respectively. Also given are results for a 4D model which does not in- 
clude water molecules 1591 . The experimental assignment is shown as well (Vnht from Ref. 1591 . 
the other modes from Refs. (7l[8l|9][l4]|34]|35l|36l|37l[38l). 



mode 


4D (AT) 


4D 6D 


exp. 


VnH2 


3330 


3297 3280 3215 


VC2=02 


1758 


1796 1792 


1716 


5nH2 


1719 


1718 1708 


1665 


VC4=04 


1645 


1709 1702 


1665 



(i) to the bending overtone 25nh2 which is by far the strongest coupUng, (ii) to the 
combination tone between the 5nh2 ™d the Vc4=04 modes. In Table|5]we also give 
the coupUngs for the isolated AT case. Notice that here only the 25nh2 overtone is 
strongly coupled. In other words, the presence of water establishes as new coupling 
channel. And, going back to Fig.|2]it requires a water molecule at the C4=04 site 
which mixes the 5nh, and the Vc4=04 modes. 

The fundamental transition frequencies obtained from this 4D anharmonic force 
field are given in Table 6 which besides the experimental values contains frequencies 
calculated for a 6D model which additionally includes most strongly coupled water 
stretching and bending modes at the N6-H site 1591 . Inspecting 4D and 6D cases we 
notice that the effect of explicit inclusion of water modes is only modest especially 
in comparison to the isolated case (4D(AT), see also Table [T]). Given the simplicity 
of the model, the agreement between theory and experiment is rather reasonable 
with deviations being about 2% except for the Vc2=02mode whose frequency is 
about 4% below the experimental value. Perhaps this is not very surprising as the 
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C2=02 mode is close to the thymine Nl site where in DNA the base is hnked to the 
backbone. 



3 Experimental Section 
3.1 Methods 

A-T DNA double strand oligomers with sodium counterions and a length of 20 base 
pairs were obtained from Biotherm, and were dissolved in water and dried on a CaF2 
window at 293 K in an atmosphere of 52% relative humidity (saturated solution of 
NaHS04.H20 at 20° Celsius |[60l). This results in DNA samples with approximately 

4 to 6 water molecules per base pair iJTll (sample thickness ~ 6.5 /im). It has been 
reported that under these conditions AT DNA oligomers adopt the B'-form |[35ll . 
Femtosecond time-resolved IR pump-probe experiments were performed with two 
independently tunable femtosecond pulses generated by parametric conversion pro- 
cesses pumped by a regenerative Ti: sapphire laser system (800 nm, repetition rate 
1 kHz, pulse duration 100 fs) (611. The central frequency of the pump pulse was 
varied from 1630 to 1760 cm^'and the probe was centred around 1650 cm^'or 
3200 cm^' . The cross correlation between pump and probe pulses had a temporal 
width of 130 fs (FWHM). With the used pump pulse energy of 1 /xJ approximately 
2% of the AT base pairs in the sample volume were excited. After interaction with 
the sample, the probe pulses were spectrally dispersed and detected with a HgCdTe 
detector array (resolution 5 cm^' ). 



3.2 Experimental Results 

The absorption between 3050 and 3600 cm^' (see Fig.|4^ and|5] solid line) is domi- 
nated for more than 85% by the broad OH stretching absorption of water molecules. 
Therefore, it is not possible to directly determine the Vnh, stretch absorption fre- 
quency of the DNA oligomer from this absorption spectrum. In the fingerprint re- 
gion, the absorption of Vc2=02 is located at 1716 cm^' and the combined absorption 
of 5nh2 and Vc4=04 peaks at 1665 cm"' (see Fig.|li, solid fine) EHIHEIlESlil- 
Both the Vc4=o4and the 5nht vibration absorb at 1665 cm^'and therefore can- 
not be excited separately in our experiment. Figure |5] shows results of femtosec- 
ond pump-probe experiments with excitation in the fingerprint region and prob- 
ing between 3050 and 3600 cm^' . Excitation with a broad pump pulse at 1740 
cm^'(FWHM 170 cm^') leads to an instantaneous spectrally narrow response 
around 3200 cm^' . Furthermore, a spectrally broad response over the entire range 
of 3050 cm^' to 3600 cm^' is seen to increase on a picosecond time scale (see Fig. 
Is). At 1740 cm^4he pump pulse mainly excites the Vc2=02 stretching vibration. 
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3600 3200 2800 1800 1700 1600 
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0.4 ps^/ ' 
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1 .4- pS s mmm^ 


pump 












ground state 



Fig. 4 (a) Absorption spectra of AT DNA oligomers prepared in 52% relative humidity (black solid 
line), neat water (dashed line), and AT DNA oligomers dried for two days in a N2 atmosphere (gray 
solid line, (b) Scheme describing the two color IR pump-probe detection of the Vnh, stretching 
vibration in AT base pairs. Because of the anharmonic coupling between the Vc4=o4 . ^nh, , and 
Vnht vibrations, the Vnht =0^ 1 transition is bleached upon excitation of the Vc4=04and Snht 
modes. Polupation of the excited state levels will locally heat the molecule, inducing a shift of the 
hot ground-state Vnht transition. 



Given the photon energy and pulse intensity, two-photon excitation of vibrations 
around 3300 cm^' is unlikely by the pump pulse. As a consequence, both the instan- 
taneous and the increasing broad negative signal must result from exciting vibrations 
in the fingerprint region. The broad negative signal, which becomes positive above 
3530 cm^ ^ , is known to correspond to the OH stretching vibration of hot bulk water. 
Excess energy in low frequency vibrations of water (e.g. librations) weakens the HB 
strength, resulting in an increase of the OH stretching force constant, and therefore 
a higher OH stretching frequency 12] |62l |63l . The instantaneous narrow response 
around 3200 cm^', after excitation at 1740 cm^', decreases in time and should 
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3600 3200 

Wavenumber (cm"') 



2800 



Fig. 5 Absorption spectrum of AT DNA oligomers around 3300 cm^' (solid line) and absorbance 
difference spectra for several pump-probe delay times after excitation at 1740 cm^' (FWHM 170 
cm^' ). The picosecond OH stretching response of water ranges from 3600 to 3050 cm^' . The 
spectrum at ps delay time were obtained by averaging from -200 to 200 fs to eliminate nonab- 
sorbing signal contributions. 




Fig. 6 (a) Contour plot of the absorbance change as a function of delay time on a log scale after 
excitation at 1630 cm^' (FWHM 160 cm^' ). On the picosecond time scale the water response is 
visible, (b) Contour plot of the absorbance change after subtraction of the monoexponential 13 ps 
rise time of the hot water response. The negative peak around 3215 cm^' is clearly visible for early 
delay times. 



therefore correspond to a different process. Figure |6] shows absorbance changes in 
the range between 3050 and 3250 cm^', upon excitation at 1630 cm^'(FWHM 
160 cm^' ), before and after subtraction of the spectrally broad 13 ps component of 
the hot water formation, obtained from a global fit. At 1630 cm^' the pump pulse 
mainly excites the 5nht and Vc4=o4 vibrations. An instantaneous bleach signal is 
observed with a maximum at 3215 cm^^ and a width of 50 cm^' , which decays on 
a subpicosecond time scale. The perturbed free induction decay of this band gives a 
total dephasing time T2 of 0.5 ±0.1 ps, which corresponds to a homogeneous line 
width of 21 ±5 cm^' . This indicates that the origin of the observed 50 cm^' width 
of the bleaching band is not caused by a single homogenously broadened absorption 
line. The maximum of the instantaneous response at 3215 cm^' decays with a 0.6 ± 
0.2 ps time constant. Around 3130 cm^' a positive signature seems to be present for 
early delay times evolving into a negative band with a rise time of about 0.4 ps, that 
decays with a time constant of 1.4 ± 0.4 ps. The time constants characterizing the 
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kinetics for various pump-probe wavelength combinations are summarized in Table 

m 

In order to identify the origin of the instantaneous bleaching signal at 3215 
cm^' we compared transients at 3215 cm^' after excitation at 1630 cm^^ and 1730 
cm^' , respectively, in Fig. [7] Excitation at 1730 cm^^ results in a w 7 times weaker 
and barely visible instantaneous negative signal around 3215 cm^' . In contrast, ex- 
citation around 1630 cm^' result in a pronounced instantaneous signal around 3215 
cm^' . Comparing the two transients at 3215 cm^' , we observe that, excitation at 
1630 cm^' leads to a three-exponential decay of the bleach with 0.6 ps, 3 ps, and 
13 ps. A transient with similar time constants (0.9 ps, 4 ps, and 13 ps), but very 
different amplitudes is observed after pumping at 1730 cm^' , where mainly the 
Vc2=02 stretching vibration is excited. Thus, excitation in the spectral range of the 
Vc4=04/ 5nh7 vibrations results in significant instantaneous signals around 3215 
cm^' . 

Experiments where both the pumping and probing takes place in the fingerprint 
region are presented in Fig.|l] In Fig.|8^ the AT DNA oligomer sample was excited 
at 1 760 cm^ ' and probed between 1 605 and 1 740 cm^ ' . The pump-probe spectrum 
shows a negative band at 1725 cm^' and a positive band at 1685 cm^' . The band 
at 1725 cm^' decays mono-exponentially with 0.9 ±0.1 ps, while the band at 1685 
cm^' decays with 0.7 ± 0. 1 ps. Transients at these frequency positions are presented 
together with their simulations in Fig. |9] The positive signal can be assigned to 
the Vc2=02 =1^2 transition. A similar lifetime was obtained by Zanni et al. for 
measurements on G-C DNA oligomers with excitation in the same frequency region 
II20II . The difference between the 0.7 ps excited state lifetime and the 0.9 ps ground 
state recovery time signals is that the Vc2=02 energy is first converted into excitation 
of lower frequency modes. As a consequence of this, the ground state absorption 
frequency is shifted due to anharmonic coupling to these lower frequency modes, 
and will not recover before these modes loose their excitation energy. 

Results for excitation of the AT DNA oligomer at 1630 cm^', presented in 
Fig. [8J3, show bleaching signals at the ground state absorption positions of the 
Vc2=02 vibration (1716 cm^'), the Vc4=o4aiid 5nh2 vibrations (both 1665 cm^'). 



Table 7 Time constants of transients, r: rising signals; d: decaying signals. 



mode 


pump (cm ' ) (FWHM) probe (cm 


' ) time constants (ps) 


VNH2 


1630 (160) 


3215 


(d) 0.6 ± 0.2 // (d) 3.0 ± 1.5 // (r) 13 ± 2 


VH20" 


1630 (160) 


3130 


(r) 0.4 ± 0.2 // (d) 1.4 ± 0.4 // (r) 13 ± 2 


VnHt 


1730 (90) 


3215 


(d) 0.9 ± 0.4 // (r) 4.0 ± 1.5 // (r) 13 


VC2=02 


1760 (100) 


1725 


(d) 0.9 ±0.1 


VC2=02 


1760 (100) 


1685 


(d) 0.7 ±0.1 


VC2=02 


1630 (130) 


1720 


(d) 2.4 ± 0.2 


VC4=04/ SnHj 


1630 (130) 


1665 


(d) 0.4 ±0.1 //(d) 1.4 ±0.4 




1630 (130) 


1650 


(d) 0.2 ±0.1 //(d) 1.0 ±0.2 




1630 (130) 


1640 


(r) 0.6 ±0.1 


VC4=04/ 5nHt 


1630 (130) 


1625 


(d) 0.5 ±0.1 



° Suggested assignment (see text) 
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Fig. 7 (a) Contour plot of the absorbance change as a function of delay time after excitation at 
1630 cm^' (FWHM 160 cm^' ) from 3140 to 3250 cm^' . The instantaneous negative signal at 
3215 cm^' as well as the dynamics of the hot water vibrations on a picosecond time scale is clearly 
visible, (b) Contour plot of the absorbance change after excitation at 1730 cm^' (FWHM 90 cm^' ). 
The signals are much weaker. The picosecond dynamics of hot water vibrations is clearly visible, 
but the instantaneous response around 3215 cm^ ' is negligible. 




Fig. 8 (a) Contour plot of the absorbance change as a function of delay time after excitation at 1760 
cm^' (FWHM 100 cm^' ). (b) Contour plot in the same frequency range (1605 to 1740 cm^' ) after 
excitation at 1630 cm"' (FWHM 160 cm"' ). 



and the water Sn^o vibration (1650 cm"' ) ||2l[3Tl|22]|63]. The Vc2=02 bleach recov- 
ers exponentially with 2.4 ± 0.2 ps on a significantly longer time scale than after 
direct Vc2=02 excitation. This shows that there are at least two different energy re- 
laxation pathways in DNA involving the Vc2=02 vibration. For the Vc4=o4/ ^nht 
vibrations biexponential recoveries were observed with 0.4 ± 0. 1 ps and 1 .4 ± 0.4 
ps, and for the water 5h,o with 0.2 ± 0. 1 ps and 1 .0 ± 0.2 ps. For the bending vibra- 
tion ShjO of water molecules in bulk water a lifetime of 170 ± 30 fs 1631 has been 
reported, which agrees with the fast component observed here for ShjO ■ Instanta- 
neous increased absorption signals below 1640 cm"' , decay with a time constant of 
0.5 ± 0. 1 ps, matching the fast component of the Vc4=o4/ bleach recovery. In 



16 



Ktihn et al. 




Fig. 9 Transients at 1685 cm '(circles) and 1725 cm '(squares) after excitation at 1760 
cm"'(FWHM 100 cm"'). Solid lines represent fits witli 0.7 ± 0.1 ps and 0.9 ± 0.1 ps time 
constants for 1685 and 1725 cm"' , respectively. System response (dashed lines). 



addition, a positive band appears around 1640 cm" , rising in 0.6 ±0.1 ps. Signals 
at this spectral position have been assigned to the bending vibration of hot water 
molecules ll63l . 

The dynamics of experiments with different pump frequencies in the fingerprint 
region, are compared to confirm that the bleach band at 3215 cm"' originates from 
the Vnh2 stretching vibration of adenine and not from water. First, we point out, 
that even the fastest recovery time of the bleach signal at 3215 cm"' after excitation 
in the fingerprint region (where the water 5h,o also absorbs), is three times slower 
than the reported lifetime of the water ShjO vibration 16311 . Second, the width of 
the bleaching band at 3215 cm"' of 50 cm"' is considerably narrower than the 200 
cm"' estimate for the 5htO overtone band of bulk water 1641 . Having assigned the 
3215 cm" ' bleaching signal to the AT DNA oligomer Vnh2 vibration, we now com- 
pare the dynamics at this frequency for excitation at 1630 cm"' vibration and 1730 
cm" ' , which correspond to the absorption bands of both the Vc4=04 and 5nh2 vi- 
brations, and the Vc2=02 vibration, respectively. 

For both excitation frequencies the picosecond dynamics, shown in Fig.|7] can be 
modelled by the same 13 ps time constant, corresponding to the rise of hot water sig- 
nal due to energy transfer from the DNA oligomer to the water molecules. However, 
the initial sub-picosecond dynamics are markedly different for these two excitation 
frequencies, further confirming that both signals can not be due to the overtone 
excitation of the water Snht ■ Although the absorption at both frequencies is compa- 
rable, the signal strength after excitation at 1630 cm" ' is in fact seven times stronger 
than after excitation at 1730 cm"' , which indicates that the Vc2=02 vibration cou- 
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pies substantially weaker to the Vnh2 vibration than both the Vc4=04 ™d the 5nh2 
vibration. 



4 Discussion 



A model summarizing the theoretical and experimental results is given in Fig. |4t). 
The calculated couplings of the Vc4=o4 - 5nh2 - ^nd Vc2=02 to the Vnht vibration in- 
dicate that excitation of either of these modes should result in a bleaching signal due 
to the shifting of the Vnht =0^1 transition. The bleach signals in Fig. [7] around 
3215 cm^' agree with this theoretical result. Furthermore, the force constants in Ta- 
ble |5] predict that excitation of the Vc4=04 and 5nht vibrations should result in a 
larger shift of the Vnht vibration than of the Vc2=02 vibration. This is confirmed by 
the data in Fig.|7] Experimentally, one cannot distinguish between the contributions 
of the Vc4=o4aiid modes. From the force constants in Table |5] however, we 
conclude that excitation of the 5nht vibration is expected to have the most sub- 
stantial effect on the shift of the Vnht vibration. The agreement between theoretical 
predictions and IR pump-probe measurements allows us to assign the bleaching 
signal at 3215 cm^'at least partially to the symmetric Vnht vibration of adenine. 
This absorption band lies about 100 cm^' lower in energy than the same mode in 
modified adenine-uracil Watson-Crick base pairs in solution |2Ti. The lower fre- 
quency of the hydrogen-bonded Vnh, vibration in DNA films compared to single 
AU base pairs in CDCI3 solution can be rationalized by significant interactions with 
neighbouring base pairs and water molecules, that weaken the force constant. The 
theoretical results further show, that inclusion of water molecules leads to a HB be- 
tween a water molecule and the NH2 group of adenine, and therefore a coupling of 
the adenine Vnht vibration and of the water bending vibration. Due to this coupling, 
this water molecule could act as a primary energy sink in energy disposal by DNA. 
Since the amount of water molecules and the coupling of water molecules to the 
nucleic acids is different for the major groove and the minor groove, one would ex- 
pect different energy relaxation pathways with deviating time constants for energy 
flow from water to DNA and DNA to water, for both sides. Energy redistribution 
processes in DNA itself provide complex relaxation patterns as presented for the 
Vc2=02 vibration. 

In summary, the presented results demonstrate the capacity of combining IR- 
pump-probe methods with calculations on microsolvated base pairs to reveal infor- 
mation on and identify hidden vibrational absorption bands. The simulation of real 
condensed phase dynamics of HBs, however, requires to take into account all intra- 
and intermolecular interactions mentioned in the Introduction. As far as DNA is 
concerned, Cho and coworkers have given an impressive account on the dynamics 
of the CO fingerprint modes ll22l l23l l24l l25l . Promising results for a single AU pair 
in deuterochloroform 11211 have been reported recently using a Car-Parrinello based 
QM/MM scheme IMl- 
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